& , 



RAMSCRAM — A Flexible 
Ramjet/Scramjet Engine 
Simulation Program 


t 

) 


Leo A. Burkardt and Leo C. Franciscus 
Lewis Research Center 
Cleveland , Ohio 


Prepared for the 

s^^Tv^a 085 Turbi " e Aeroengine Congress and Exposition 
sponsored by die American Society for Mechanical Engineers 
Brussels, Belgium, June 11-14, 1990 * 


NASA 


(NAS*-TM-10?4®> 1 ) KAMSC°AM: a H_ c X I r i l ! 
HAMJET/SCPAMjr T cNG I N r MMUtATlON P o 00RAM 
(NASA) ti p CSCL ?1F 


N 1 ? 0-14^35 



0 3/07 


Unc 1 as 
0?5445V 




1 


! 1 

r 

i I 


RAMSCRAM - A FLEXIBLE RAMJET /SCRAMJET ENGINE SIMULATION PROGRAM 


Leo A. Burkardt and Leo C. Franciscus 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

With the resurgence of interest in high super- 
sonic and hypersonic flight there is a need to simu- 
late airbreathing engines which may be used in this 
flight regime. To meet this requirement the RAMSCRAM 
code was developed. The code calculates one-dimensional 
flow properties at each component Interface and the 
overall performance of the engine. It uses equilibrium 
thermodynamics which accounts for dissociation and 
allows for any fuel or combination of fuels. The pro- 
gram can simulate ramjet, scramjet, rocket, and ducted 
rocket engines. 

INTRODUCTION 

In the past few years, there has been a resurgence 
of Interest in high speed flight. Such programs as the 
High Speed Civil Transport and the National Aerospace 
Plane have highlighted this Interest. Ramjet and scram- 
jet engines perform very favorably in the high super- 
sonic and hypersonic flight regimes. Ramjets can be 
effectively used from about Mach 2 to Mach 6 and scram- 
jets give reasonable performance at Mach 6 and above. 

The RAMSCRAM program described herein was deve- 
loped at NASA Lewis Research Center to meet the need 
to calculate the performance of ramjet and scramjet 
engines. The goal In developing this program was to 
obtain high fidelity simulation while maintaining ease 
of use. versatility and reasonable calculation times. 

In pursuit of this goal, performance Is calculated by 
marching from front to back through the engine flow 
path In a stepwise fashion. One-dimensional calcula- 
tions are performed at each step with efficiency fac- 
tors to account for flow path losses. Program Input 
Is made easy through the use of Fortran Namelist Input. 

Equilibrium thermodynamics is Included to account 
for high temperature real gas effects. The general 
nature of the thermodynamics subroutine, together with 
the flexible mixer and combustor components, gives the 
RAMSCRAM code the ability to simulate a large variety 
of engine/propellant combinations. Some examples 
Include: any chemical rocket or ducted rocket, multi- 

fuel ramjet engines, oxygen enriched scramjet engines. 


This paper presents a general description of the 
calculation methods used in the code and a demonstra- 
tion of some of Its capabilities through the use of 
examples . 

FLOW PATH DESCRIPTION 

Figure 1 is a schematic of the simplest flow path 
jsed in the RAMSCRAM program. Calculations are per- 
formed at each flow station shown in the schematic. 
Station 0.0 Is the beginning of the forebody (external 
compression surface) and has free stream conditions 
associated with it. Station 1.0 defines the beginning 
of the Inlet (internal compression portion). Station 
2.0 defines the throat of the inlet or the end of the 
mixer if a mixer is to be included. Station 3.0 Is the 
diffuser exit/combustor entrance location. Station 4.0 
Is the combustor exit. Station 4.5 Is the nozzle 
throat. Station 5.0 Is the nozzle exit. 

THERMOOYNAMIC ROUTINE 

Due to the high temperatures that can be reached 
in ramjet and scramjet engines, we felt that equilib- 
rium thermodynamics would be needed to maintain good 
accuracy. At the time we were developing this code, 
the CET86 (1) general equilibrium code was being modi- 
fied to run as a subroutine to the NNFP (2) turbine 
engine simulation code. It was convenient therefore 
to use this version of the CET86 code for the RAMSCRAM 
program also. 
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FIG. 1 - RAMJET STATION LOCATIONS. 
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In addition to the modifications mad* to CET86 for 
us# with NNEP, further modifications were made to allow 
more than a single fuel/oxidizer combination to be 
specified. Up to 24 reactants can be specified as pos- 
sible fuels or oxidizers during a single run with the 
appropriate combination of those reactants specified 
at each location In the engine. This permits calcula- 
tion of propulsion systems such as ducted solid rocket. 

Forebod v 

A conical , two-dimensional wedge or no forebody 
can be specified. Figure 2 presents a schematic of the 
forebody. 

In the case of a cone, the conical haif angle Is 
specified and a conical shock Is calculated. The coni- 
cal flow field Is Integrated across the station 1.0 
Inlet height and a mass weighted average Is calculated 
to obtain one-dimensional flow conditions at that 
location. 

If a two-dimensional wedge Is specified, up to 20 
ramps can also be Included. This allows an Isentropic 
compression surface to be approximated. The shock off 
each ramp Is calculated. If the combination of Mach 
number and ramp angle at any ramp will not support an 
attached oblique shock, then a terminal normal shock Is 
calculated at that point and the subsonic flow field 
Is assumed to persist down the rest of the forebody. 

The Inlet can be specified to be either started or 
unstarted. If it is started, the conditions after the 
oblique shocks specify the flow field at station 1.0. 

If the inlet Is specified as unstarted, then a further 
normal shock calculation is performed as part of the 
forebody calculation to give the properties of the flow 
a station 1 .0. 

Inlet 

In the simplest case, the Internal inlet consists 
of everything from station 1.0 to 3.0. The flow char- 
acteristics at each station are based on those calcu- 
lated at the previous station and the flow area and 
efficiency terms given for the current station. The 
efficiency terms are always relative to conditions at 
the previous station. If no efficiency factor is pro- 
vided an Isentropic calculation is performed. 

The flow field at station 1.0 defines what is com- 
ing into the Inlet. An efficiency factor can be applied 
to the flow at this point to account for any losses not 
accounted for In the forebody calculations. The effi- 
ciency factor can be In one of three forms: total pres- 
sure recovery, kinetic energy efficiency or process 
efficiency. Process efficiency Is similar to kinetic 
energy efficiency, except that it Is referenced to dif- 
fuser exit static conditions rather than total condi- 
tions. If desired, it Is entirely possible to leave 
the forebody undefined and account for ail losses at 



this point using the efficiency factor. Two of three 
parameters, airflow, Mach number or flow area, are 
Input to permit calculation of th# flow conditions and 
geometry at this location. If only airflow or flow 
area are given, then the Mach number 1$ assumed to be 
th# terminal Mach number from the forebody calculation 
or the free stream Mach number If no forebody Is 
Included. If flow area and airflow are Input, then the 
subsonic or supersonic solution Is chosen to match the 
condition of the air Just upstream of this station. It 
can also be forced to either solution by entering a 
subsonic or supersonic Mach number which will be used 
as a first guess In the calculation procedure. 

Station 2.0 can be used to determine the flow 
field at an Intermediate point in the Inlet, such as 
the throat of a ramjet inlet. The same efficiency fac- 
tor and either Mach number or flow area are input to 
calculate the flow conditions. If desired, a normal 
shock can be included in the calculation at this point 
when the inlet is started. 

The flow conditions at station 3.0 are calculated 
In the same manner as descrloed for station 2.0, witn 
the exception that there is no provision for a normal 
shock calculation. 

Mi xer 

The inlet can be configured to perform as a mixer 
in the case of a ducted rocket or ejector ramjet simu- 
lation. In this case, station 1.0 becomes the air- 
stream entrance to the mixer and station 2.0 becomes 
the mixed plane. If it is not convenient to use sta- 
tion 1.0 as the mixer entrance, then another station 
can be Included to follow station 1.0. This station 
is designated station 1.5 and it has the same charac- 
teristics as described for station 2.0 in the inlet 
section above. The rocket stream enters the mixer at 
station 1.2. 

The mixer model is defined by a control volume as 
shown In Fig. 3. Streams A (air stream) and 8 (rocket 
exhaust) enter the control volume at flow station 1.5 
and 1.2, respectively, and mix Inside the control vol- 
ume The flow is assumed to be fully mixed and departs 
the control volume at station 2.0. 

The mixed flow conditions at station 2.0 are cal- 
culated by simultaneously solving the equation of state 
and the conservation of mass, momentum and energy equa- 
tions. This leads to a quadratic equation. The two 
solutions represent a low and a high entropy gain. The 
high entropy solution always results in subsonic flow 
and the low entropy solution normally results In super- 
sonic flow. The code allows the choice of etther solu- 
tion. Normally the subsonic solution would be chosen. 



FIG. ) - MIXER MODEL . 
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The mixer calculation allows the Inclusion of 
efficiency factors to account for bosses In the mixing 
section. These factors Include- a momentum loss term 
for either Incoming stream, a wall friction loss term 
and a total or static pressure loss term. The pressure 
loss term is referenced to mixed pressure before the 
loss Is calculated. 

The mixing section does not have to be constant 
area. If a diverging or converging area Is defined, 
then a pressure-area term must be accounted for on the 
walls. This term Is approximated by assuming that the 
pressure acting on the wall Is the average of stream A 
and mixed stream static pressures. 

The use of equilibrium thermodynamics also accounts 
for any burning occurring In the mixing process by main- 
taining equilibrium chemistry at the mixed plane. 

Combustor 

The combustor model Is essentially the same as the 
mixer model. In this case, the fuel stream replaces 
the rocket exhaust stream of the mixer. Two fuel 
streams are permitted to allow the Injection of two 
different fuels simultaneously. Also, up to ten 
sequential fuel Injection stations are permitted along 
the length of the combustor. A ramjet combustor Is 
represented by specifying the subsonic solution and a 
scramjet combustor is represented by specifying the 
supersoni" solution. Using this model for the combus- 
tor allows the momentum of the fuel to be accounted for 
as well is its energy. The fuel momentum can be a sig- 
nificant contributor to scramjet thrust. 

Incomplete combustion is represented by Including 
••Inert fuel" in the fuel composition. Currently "Inert 
hydrogen" and "Inert methane" are available, but others 
can be easily added as needed. Other combustor effi- 
ciency factors are similar to those used for the mixer. 

Nozzle 

The nozzle calculation Is an isentropic expansion. 

A throat calculation is provided whenever the velocity 
coming out of the combustor Is subsonic. The nozzle 
can be allowed to expand fully to free stream pressure 
or It can be terminate at a given area. No distinction 
Is made between Internal and external expansion sections 
of the nozzle. Two efficiency factors are provided, a 
velocity coefficient and a thrust coefficient. 

EXAMPLE SIMULATIONS 

Ducted Rocket 

A ducted rocket was simulated and an analysis of 
its performance versus air flow to propellant flow 
ratio was completed. The rocket burns a composite 
solid propellant mixture of ammonium perchlorate oxi- 
dizer and ethylene oxide fuel. The net thrust was 
held constant while the air to propellant bypass ratio 
was varied. As the bypass ratio was varied, the mix- 
ture ratio of the propellant was also varied to give 
enough excess fuel to provide a stoichiometric fuel to 
air ratio. The specific Impulse and air capture area 
of the rocket versus bypass ratio Is plotted In Fig. 4. 
The capture area Is normalized with thrust. 

A bypass ratio of 6 was arbitrarily chosen as the 
maximum simulated. At some point, the solid propellant 
mixture ratio becomes such that It cannot support com- 
bustion. The specific Impulse Increases linearly with 
bypass ratio. However, as shown, the air capture area 
also Increases with bypass ratio. This offsets the 
gain In performance, as the Inlet and combustion duct 
will Increase In size and weight. The ultimate compro- 
mise between the two offsetting trends must be accom- 
plished by doing a vehicle analysis using the Informa- 
tion provided by the engine performance code. 
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FI6. 4 - DUCTED ROCKET PERFORMANCE. 




FI 6. 5 - PERFORMANCE COMPARISON OF CANDIDATE RAMJET FUELS. 


Ramjet Fuel Comparison 

A ramjet was operated at Mach 5 using three dif- 
ferent fuels. The general thermodynamics routines used 
in this program make this type of comparison a very 
simple procedure, requiring only a change In the speci- 
fied fuel to air ratio and the fuel type. Figure 5 
illustrates the performance of the engine with each 
fuel . 

Figure 5(a) presents Icn and Fig. 5(b) presents 
volumetric Icn, whlcn Is defined as thrust divided by 
volumetric flow rate. The three fuels represented are 
JP 4 , liquid methane and liquid hydrogen. As shown, 
hydrogen gives by far the best performance based on 
mass flow rate, but it also gives by far the worst per- 
formance based on volumetric flow rate because of its 
very low density. The structural weight of the vehicle 
required to accommodate the fuel and the larger drags 
associated with the greater volume of hydrogen or meth- 
ane play an important part In the tradeoff between 
fuels. Here again, the performance provided by the 
engine code would be used In a vehicle analysis to 
determine the best fuel for specific mission 
requirements . 

Scramjet Sensitivity Analysis 

A sensitivity analysis was performed for a 
Mach 10 scramjet. The Inlet kinetic energy efficiency 
(n<£), combustion efficiency (nfl)* an d nozzle thrust 
coefficient (Cpq> *e r e varied. The baseline values 
used were: 0.97 for n<£. 0.95 for 09 , and 0.96 for 

Cfg- Figure 6 presents the results of this analysis. 

Each parameter was varied over a range of from 
approximately 0.9 to 1. Kinetic energy efficiency was 
varied with two different downstream constraints. In 
one case, the diffuser exit/combustor Inlet geometry 
was held constant. In the other case the diffuser 
exit/combustor inlet Mach number was held constant. In 
the first case hke had only a small effect on engine 
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CONCLUDING REMARKS 



EFFICIENCY PARAMETER 
FIG. 6 - SCRAMJET SENSITIVITY ANALYSIS AT MACH 10. 


performance because there was a self compensating 
effect of the burner Mach number decreasing as hke 
decreased. This resulted In lower Rayleigh losses in 
the burner. Combustion efficiency also had only a 
small effect on the performance over its range of vari- 
ation. Both n<E with constant burner Mach number and 
CpQ had large effects on the engine performance. This 
is because they both directly affect the gross thrust 
on about a one-for-one basis. Net thrust is only about 
15 percent of the gross thrust and therefore any varia- 
tion in gross thrust Is magnified in the net thrust and 
specific impulse values. 


Thtrt has been a resurgence of Interest In hyper- 
sonic propulsion. Air breathing flight In this regime 
requires the use of ramjet and scramjet propulsion sys- 
tems. Vehicle and engine performance and tradeoff 
studies require a reasonably accurate, versatile and 
fast running propulsion simulation program. 

The RAHSCRAM code fulfills this requirement. The 
general form of the geometry Inputs allows the simu- 
lation of many engine configurations. Component 
efficiency Inputs make It easy to study component per- 
formance effects on overall engine performance. And 
the general nature of the equilibrium thermodynamics 
routines allow the specification of nearly all conceiv- 
able propellant combinations. 
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